Comment on "Time-Dependent Density-Matrix 
Renormalization Group: A Systematic Method 
for the Study of Quantum Many-Body Out-of- 
Equilibrium Systems" 

In a recent Letter Q, Cazalilla and Marston (CM) 
proposed a time-dependent density-matrix renormaliza- 
tion group (TdDMRG) algorithm for accurate evaluation 
of out-of-equilibrium properties of quantum many-body 
systems. An interesting result they obtained is that for a 
point contact junction between two Luttinger liquids, a 
current oscillation develops after initial transient in the 
insulating regime (Fig. 2 in Ref. 1). They attributed 
this oscillation to a non-linear response of the Luttinger 
liquid. Here we would like to point out that (a) the os- 
cillation they observed is an artifact of the method; (b) 
the TdDMRG can be significantly improved by extend- 
ing the definition of the density matrix to adapt the non- 
equilibrium evolution of the ground state. 

In the TdDMRG scheme of CM, it is the equilibrium 
ground state that is targeted and used in the construc- 
tion of the reduced density matrix. At the beginning of 
the evolution, the state does not deviate much from the 
ground state \^o), and results obtained with that ap- 
proach are accurate. However, the time evolution of the 
ground state evokes excited states in a non-equilibrium 
system, so the long time behavior of the state becomes 
very poor and could not be substantially improved by 
keeping more states. This can be clearly seen from Fig. 
1, where the numerical results obtained with that ap- 
proach (i.e. the Nt — curves) for the two models in 
Ref. [j] begin to deviate from the exact ones at t ~ 25. 
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FIG. 1: Current for (a) a quantum dot and (b) a junction with 
V = 0, defined by Eqs. (2) and (5) in Ref. |l], respectively. 
The parameters are the same as in Ref. ^. 

To retain the information on the relevant excitation 
states, we propose to define the density matrix from the 
time dependent wavefunction as 

N t 

P = J2^(U))(^(U)l 5>^ = 1, (1) 

i=0 i 

where | 'J' (to)) = h^o), is the starting time at which a 
bias is imposed and Nt is the number of equal time inter- 
vals within the whole evolution time to be evaluated. The 
reduced density matrix for the system is defined by trac- 
ing over the degrees of freedom of the environment. Both 


the time evolution of the wavefunction and the newly de- 
fined reduced density matrix should be evaluated from 
the beginning of the DMRG truncation, not just for the 
final lattice system. The density matrix used by CM cor- 
responds to the Nt = case. In our calculation, we took 
on = 1/2 for i = and l/2N t otherwise. Empirically, we 
find that Nt > 6 is generally needed. Fig. 1 shows our 
results with Nt — 24. By comparison with the exact re- 
sults, we find that this new scheme significantly improves 
the large time scale behavior of the current. 
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FIG. 2: Tunnelling current for the model defined by Eq. (5) 
in Ref. [1] with V = l.lw and Sfj, = 6.25 x 10" 2 uj. 

Now let us turn to the tunnelling junction between two 
Luttinger liquids in the insulating regime. For N t = 0, 
our calculation confirms the existence of a current oscilla- 
tion around t ~ 20. However, we find that the oscillation 
depends on the number of states retained. By further ap- 
plying the new TdDMRG scheme to this system, we find 
that the oscillation does not exist at all. In Fig. 2, the 
numerical results with N t = and Nt — 24 are shown for 
detailed comparisons. Therefore, the current oscillation 
observed by CM is in fact an ar tifact of the TdDMRG [Q, 
rather than a nonlinear response of the system to the ap- 
plied bias. This also shows that in the study of physical 
properties of a quantum state out of equilibrium, it is 
important to include the relevant excitation states in the 
definition of the reduced density matrix. 
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